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The shape function describes (non-pertnrbative) Fermi motion eflects in semi-inclusive heavy flavour decay. Its 
renormalization properties are substantially dependent on the kind of ultraviolet regulator used. For example, 
the double logarithm that appears at one loop is larger by a factor 2 in dimensional regularization than in lattice 
regularization. We show that factorization of long-distance effects inside the shape function is achieved with any 
regulator considered. 



1. Introduction 

In general, the shape function |l| is introduced 
to describe the decay of a hadron Hq containing 
a heavy flavour Q into an inclusive hadronic state 
X with a large energy and a small invariant mass, 
plus non-QCD partons, i.e. 

Hq -^ X+non - QCD partons {mx < Ex) .(1) 

Specifically, the selected kinematics is 



-j^-O (Aqcd) , 



Ex :$> Aqcd- 



(2) 



In a more formal language, the conditions (0) cor- 
respond to the limit 



E 



X 



cx), mx 



OO, with — — — > const. (3) 
Ex 

The limit (ra) implies the infinite mass limit of the 
heavy flavour, 



rriQ 



(4) 



as rriQ > Ex- The heavy quark can then be 
treated in the Heavy Quark Effective Theory 
(HQET) 0. The divergence oimx - even though 
it is slower than the one of Ex - implies that 
the final hadronic state can be replaced by a par- 
tonic one and that perturbation theory (PT) can 
be applied; however, the decay has also a non- 
perturbative component - to be factorized in the 
shape function - which we introduce with the fol- 
lowing model. 

We identify a hard subprocess in (|]) consisting 
in the fragmentation of the heavy flavour Q, 



Q ^ X + non — QCD partons. 



(5) 



X differs from X in that it does not contain the 
light valence quark(s) of Hq. The momentum of 
Q can be written as [0 ^ 



Pq = mnv + k', 

where mn is the mass and 

yf" = (1; 0,0,0) 



(6) 



(7) 



is the velocity of Hq , which we take at rest with- 
out any loss of generality. The momentum of the 
light degrees of freedom in Hq is —k'. The point 
fc' = is the elastic one where all the initial light 
partons have soft momenta. It is natural to as- 
sume that HI 



(8) 



k'^ ^ k'_ ^ k'^ ^ O (Aqcd) - 



The distribution of the fc'-momenta is clearly 
non-perturbative and it is related to the well- 
known Fermi motion of the heavy quark inside 
the hadron M. Let us denote by q'^ the momen- 
tum carried away by the non-QCD partons, taken 
along the +z axis; the flnal hadronic system flies 
along the minus direction. X has a mass 



m| = (g + fc')^ 



mJi+2Ex fcV 



+i-^k'_+k",{9) 



2E 



X 



where Q^ is the momentum available to the final 
partons q 



Q 



mnv 



{niH -9o; 0,0,-173) 



^ We use mj] instead of mg to define the mprnentum (H) ; 

this amounts to a shift in the range of k', [Bl . 

^Tho Hght-cone components are defined as a± = ao it 03. 
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(10) 



with v± — (1; 0, 0, ±1). The sizes of the terms in 
the third member of eq. (bl) are 



2Exk\ ^ 0{ExKqcd): 



m 



X iJ 



2Ex 



k'_ ~ k'^r^O{KlcD) 



(11) 



We now perform two different approximations: 

1. We hnearize the problem by dropping the 
k'^ term in eq. (H), i.e. we describe the 
final hadronic system with the HQET or, 
in geometrical language, with a Wilson line 
Ojff the light cone; 

2. We drop the term m\/ {2Ex) k'__ in the last 
member of eq.(^. Q^ is replaced with a 
vector lying exactly on the light cone, in 
the minus direction. Neglecting these mass 
(virtuality) effects, the final hadronic sys- 
tem is described by the LEET [^ or, in ge- 
ometrical language, by a Wilson line on the 
light cone. 



Therefore, we have 



X 



i-X 



2Exk'. 



(12) 



Equation (O) is the main result of this section 
Let us comment on it. 



X 



m^ depends on a single 



light-cone component, fc^in our reference frame 
Because of eqs. (0) and (0), the two terms on the 
r.h.s. of eq. (p^) are of the same order: thus rrvL is 
affected in a substantial way by the distribution 
of the k'_^ momenta and cannot be considered con- 
stant. 

The rate of the elementary process (g|) contains 
perturbative corrections of double-logarithmic 
kind: 



J log" [I -x] 



(0<fc<2n-l), (13) 



where ^ 



x = l- 



X 

^X 



[E- - Ex). 



(14) 



Since rn?x and w?^ are of the same order of mag- 



X 



nitude, the limit (g) implies that the threshold 
region (also called large-x region) is approached: 



1" 



(15) 



The perturbative corrections ( p^ ) — enhanced 
in the threshold region — are large and radi- 
cally modify the tree-level distribution |0]. The 
physical distribution for (Q) is obtained by con- 
voluting the perturbative corrections of the form 
(O) with the primordial fc(,_-distribution: this is 
the way non-perturbative effects enter the game. 
We conclude that the process has a substantial 
non-perturbative component related to the k'_^_- 
distribution. 

Let us also present another way to establish 
the non-perturbative component in the decay (|l|), 
which is a critical analysis of dynamics using only 
perturbation theory. We neglect all the bind- 
ing effects (confinement, Fermi motion, etc.) and 
consider an isolated on-shell heavy quark Q, i.e. 
we take fc' = Q Q decays into a massless quark 
q plus non-QCD partons: 



q + non — QCD partons. 



(16) 



We now consider the emission of a soft gluon, with 
momentum components of the order of the QCD 
scale: 

k+^k^^ k^ ^ O {Kqcd) ■ (17) 

The invariant mass of the final hadronic state is 



•-x 



= {p, + kf ~ 2Ex k+ ^ O [Ex Kqcd) ■ (18) 



An invariant mass of the order of (^, i.e. rather 
large, is generated by the emission of a very soft 
gluon. A kinematical amplification by a factor 
Ex has occurred as the kinematics goes from 

3 The plus-distribution P{x), is defined on test 

functions / {x) in the unit interval x € [0, 1] as 

jlP(x)[f{x)-f{l)]dx. 

''The physical process and the subprocess coincide in this 



time-like to light-like with the fragmentation of 
the heavy flavour. Higher-order perturbative cor- 
rections replace the bare coupling with the run- 
ning coupling evaluated at the transverse gluon 
momentum [g, 

as— >as(fc±). (19) 

For a soft gluon with the momentum ([l7|), the 
coupling is evaluated close to the Landau pole, in- 
dicating the presence of non-perturbative effects. 
We conclude again that the decay (|l]) has the non- 
perturbative component identified before. 

Let us observe that the linearization introduced 
in 1. and the light-cone limit introduced in 2. 
are not valid approximations for a hard coUinear 
gluon; by this we mean a gluon with momentum 
components of order 

k^^O{Ex), k+r^OiAgcD) 

and, in order to have fc^ r^ 0: 



O 



(V^ 



A 



QCD 



(20) 



For such a gluon, all the terms in the last mem- 
ber of eq. (g) are of the same order. Its contri- 
bution can be considered a short-distance effect, 
since the transverse momentum ( |2(]| ) is substan- 
tially larger than that of a soft gluon ( |l7|) and 
the related coupling constant is in the perturba- 
tive region. 

According to the above discussion, it is natural 
to identify the non-perturbative component in (|^) 
as the following matrix clement: 

/ (fc+) EE {Hq {v) \hl6 (fc+ - iD+) K\Hq («)), (21) 

where /i„ is a field in the HQET with velocity 
V. This function gives the probability that the 
heavy quark in the hadron has a plus virtuality 
equal to fc+, independently from the other com- 
ponents. In the following section we will critically 
analyse how the above matrix element factorizes 
the non-perturbative effects. We already found 
that the hard coUinear region (20) cannot be de- 
scribed by the shape function but can consistently 
be incorporated into a coefficient function ^ |0] . 

^The coefficient function is also called jet factor, coUinear 
factor, matching constant, hard factor and short-distance 
cross-section. 



2. The shape function in various regular- 
izations 

In order to avoid distributions and to deal only 
with ordinary functions, it is convenient to con- 
sider the light-cone function 0] 



F {k+) ^ {HQ\hlr 



1 



-h„\HQ), (22) 



' A;+ — iDj^ + ie 

from which the shape function is obtained by tak- 
ing the imaginary part, 

1 



f{k. 



-lmF{k^ 



(23) 



We want to study if the light-cone function factor- 
izes the non-perturbative effects by performing a 
perturbative computation. This involves 3 steps: 

1. To replace the hadronic light-cone function 
with a partonic light-cone function H, 



Fik+) = {Q\hl- 



c_l_ — iD^ 



le 



-K\Q), (24) 



and to perform a perturbative computation. 
We take the heavy quark Q with the mo- 
mentum (||), i.e off-shell by k'; 

2. To perform the same perturbative compu- 
tation — with the same external states — 
in the original high-energy theory, i.e. full 
QCD; 

3. To compare the results to see if the differ- 
ence is a short-distance effect or not. 

All this procedure has a meaning if we ac- 
cept the following assumption: the long-distance 
effects of perturbative kind — i.e. the lead- 
ing infrared logarithms — are able to trace 
the long-distance effects of non-perturbative na- 
ture. Therefore, if two matrix elements have the 
same perturbative long-distance contributions, 
they manifest also the same non-perturbative 
long-distance effects if computed with a non- 
perturbative technique, such as lattice QCD [|. 

®The same symbol F (k^) is used for the two different 
matrix elements, as this should not cause confusion. 
^An exception to this rule seems to be the observed differ- 
ence of the fragmentation functions of the heavy flavours, 
as extracted from electron and proton collisions |lO|. 



At the tree level, we have the expected result 

^ (25) 



F(fc+)-- = 



k'j^ — /s+ + ie ' 



I.e. 



/(fc+)*- = 6 (fc; - fc+) 



(26) 



The shape function is a spike when its argument 
matches the plus virtuality of the external state, 
fc+ — k'^, independently of the other components, 
k'_ and k'^. We now consider the double loga- 
rithm that appears at one-loop. 

Let us first discuss the case of a simple regular- 
ization cutting the space momenta: 



M<As, 



-CXD < ^0 < +00. 



(27) 



Since infrared logarithms are associated to quasi- 
real configurations, for which Zq ~ I ^ I j we expect 
this regularization to give the same double loga- 
rithm as lattice regularization (after continuation 
from Euclidean to Minkowski space). On the lat- 
tice, all the momentum components are cut off. 



< 



(28) 



where a is the lattice spacing. The result of the 
computation done with the regularization (p7|) is 



F(fc+) 



As 



1 /a 
T^ V 2 



log^ 



A, 



k + iO \ 21 ° \-k-iO 
where we defined the usual combination 



k' 



k+. 



,(29) 



(30) 



and a = asCp/T^- 

The result of the computation done with Di- 
mensional Regularization (DR) is |l^J^,|l^ Q 

' ^-?)r(i + .) 



FB{k+) = 



k + iO\ 2. 

r(l + 2e)r(l-2e) 



-k-iO 



2c 



1 



k + iO 

-log 2 



1 1, 

--e^^'H-k-^O 



2e2 



Ai 



"This result is also in 1 1 
malization procedure 1 9 



k-iO 

but we disagree with the renor 



(31) 



where e = 2 — D/2 (with D the space-time di- 
mension) and /i is the regularization scale Q. We 
find a double pole of UV nature because infrared 
singularities (soft and coUinear) are completely 
regulated by the light-quark leg being off-shell. 
The first problem we encounter is how to renor- 
malize the above expression. The point is that, to 
obtain a finite result, one has to subtract not only 
the double pole, whose coefficient is just a num- 
ber, but also the simple pole, which has a log/x/fc 
coefficient ||lj,0. It seems that the subtraction 
of the simple pole cannot be performed with a lo- 
cal counter-term, as standard textbooks on renor- 
malization teach one should do. However, if we 
blindly subtract all the poles, we obtain 



F(fc+)" 



iO 



i-a) log^ 



fJ- 



-k-iO 



.(32) 



We see that the coefficient of the double logarithm 
is 2 times larger than in eq. (ES). 

Even though we are not able to justify the 
renormalization procedure, we believe that the 
above result is correct. A hint in favour of this 
conjecture is obtained by repeating the same com- 
putation in the bare theory at finite ultraviolet 
(UV) cut-off, using a regularization similar to 
DR, i.e. cutting only the transverse momenta and 
not the longitudinal ones Fj: 



\1±\<A±, 



-CX) < A+, Z_ < oo. 



(33) 



The unbounded integration over the longitudinal 
components does not give rise to UV divergences 
and we obtain the result IHl 



Fik+)^- 



1 



k + iO 



(-a) log 



A 



-k~iO 



(34) 



The transposition of symbols should be clear: the 
renormalization point /i is replaced by the UV 
cut-off Aj_. As anticipated, there is agreement 
between the two results (|3^) and (p4|). 

The origin of the factor 2 in the difference of 
the double logarithms is easily seen by considering 



® Actually, in the last member of ( 
placed by ^■^47r exp[— 7^;] (7^ is the 



), fi should be re- 
uler constant), even 

though this rescaling does not affect DLA results. 

^''The loop measure in DR is indeed regulated as / d^l = 

1/2 J dl+dUd^-^l^. 



the integrals 

Jo ^ Jo '± 

d{As-e),e{Al-ll). 

A straightforward computation indeed gives: 



^As 



1 



■log- 



As 



Ia^ = log' 



A 



Let us analyse the domains of the 2 integrals in 
the plane {e,l±). For /as, the integration region 
is a triangle-like domain Ds limited by the curves 



li 



l± = e and e = As {Ds) ■ 



The energies and the transverse momenta in Ds 
are in the ranges (see fig.l) 



k+ < e < Ai 



k+ <l± < As 



(Ds). 



Transverse momenta do become as small as fc+ ^ 
O {Aqcd), entering the non-perturbative region. 
For /a^ , the integration region D±^ has again a 
triangle-like shape and it is limited by the curves 



li 



e and / 1 



A, 



{Da 



Assuming A_l = As , we see that D±^ contains Ds 
plus another region AZ) (see fig. 2), 

D± = DsU AD. 

In the "new" region AD 

<li < A I , 



Vfc+Aj 



A I < e < 



A2 



{AD). 



The integration over AD gives a double logarith- 
mic contribution equal to the one coming from 
the integration over Ds'. that is the origin of the 
factor 2 in the difference between /a^ and Ja^^ . 
The contribution coming from AD is a short- 
distance contribution because it always involves 
rather large transverse momenta, 



li_ > \/k+A± > A 



QCD 



(Ax » Aqcd). (35) 



Let us observe also that in AD the energy e of 
the gluon can become much larger than the UV 
cutoff, e » A^. One usually relates the cutoff 
to the hard scale Ex in the full QCD process, 
which, on the contrary, sets the upper bound for 



the gluon energies: e < Ex- It is clear that we 
are dealing with highly unphysical regularization 
scheme effects. 

We now compare with the QCD rate, from 
which a full QCD light-cone function can be de- 
fined M , given by 



F{k+)'^^'' 



1 



iO 



-f) 1-^ 



E 



X 



iO 



•(36) 



It is clear that A_l or /i is replaced by Ex. It 
is immediate to see that the double logarithms 
match for the As-regularization, while they do 
not match for the A_L-regularization. In the lat- 
ter case the "spurious" contribution from AD can 
be subtracted by means of the coefRcient func- 
tion, which includes also "true" hard coUinear ef- 
fects, finite terms, etc. While in the As case the 
coefficient function contains at most a single log- 
arithm of fc+, in the A^ case it contains also a 
double logarithm of fc+. 

3. Conclusions 

We tried to describe in a way as transparent 
as possible the problem of factorization of non- 
perturbative effects in the decay (m) by means 
of the shape function. The starting point of our 
analysis is the observation that the shape func- 
tion exhibits a different double logarithm at one 
loop, depending on the kind of UV regulariza- 
tion adopted 0. In a regularization such as DR 
— cutting only transverse loop momenta — the 
double logarithm has a coefficient 2 times larger 
than the one derived with a lattice-like regulariza- 
tion. We have explicitly shown that the difference 
between the double logarithms is related to the 
integration over a region of large transverse mo- 
menta, so it is a short-distance effect. The QCD 
rate has a double logarithm equal to the one de- 
rived in lattice-like regularization. This implies 
that the coefficient function in lattice-like regu- 
larization contains at most a single logarithm of 
fc+ (the infrared scale of the problem), while in 
DR it contains also a double logarithm of fc+ . In 
both regularizations, however, non-perturbative 
effects are factorized inside the shape function. 

In general, we observe that there are UV reg- 
ularization effects similar to those found for the 



shape function in other operators incorporating 
Wilson hnes on the hght cone |la] . 



Acknowledgements 

I would like to thank S. Catani, A. Hoang and 
G. Korchemsky for discussions. 

REFERENCES 

1. I. Bigi, M. Shifman, N. Uraltsev and A. Vain- 
shtein, Phys. Rev. Lett. 71:496 (1993); A. 
Manohar and M. Wise, Phys. Rev. D 49:1310 
(1994); M. Neubert, Phys. Rev. D 49:3392 
(1994), ibid. D 49:4623 (1994). 

2. F. Bloch and A. Nordsieck, Phys. Rev. 52:54 
(1937); L. Foldy and S. Wouthuysen, Phys 
Rev. 78:29 (1950); D. Yennie, S. Frautschi 
and H. Suura, Ann. Phys. 13:379 (1961) 
E. Eichten and F. Feinberg, Phys. Rev, 
D 23:2724 (1981); N. Isgur and M. Wise, 
Phys. Lett. B 237:527 (1990); H. Georgi 
Phys. Lett. B 240:447 (1990); F. Hussain et 
al., Phys. Lett. B 249:295 (1990); U. Aglietti 
and S. Capitani, Nucl. Phys. B 432:315-336 
(1994). 

3. U. Aglietti et al., Phys. Lett. B 432:411-420 
(1998). 

4. U. Aglietti, Phys. Lett. B 281:341 (1992). 

5. G. Aharelh et al., Nucl. Phys. B 208:365 
(1982). 

6. M. Dugan and B. Grinstein, Phys. Lett. B 
255:583 (1991); U. Aglietti, Phys. Lett. B 
292:424 (1992). 

7. J. Kodaira and L. Trentadue, SLAG preprint 
SLAC-PUB-2934 (1982) and Phys. Lett. B 
112:66 (1982), B 123:335 (1983); S. Gatani 
and L. Trentadue, Nucl. Phys. B 327:323 
(1989). 

8. D. Amati et al. Nucl. Phys. B 173:429 (1980). 

9. U. Aglietti and G. Ricciardi, Phys. Lett. 
B 466:313 (199 9); GERN-TH/2 000-071 and 
DSF-T/2000-7, |hep-ph/0003"l46| , to appear in 
Nucl. Phys. B. 

10. G. Martinelli, private communication; S. 
Frixione and M. Mangano in "Proc. Work- 
shop on Standard Model Physics (and More) 



at the LHG", Resport GERN 2000-004, p. 
254. 

11. G. Korchemsky and G. Sterman, Phys. Lett. 
B 340:96 (1994); A. Grozin and G. Ko- 
rchemsky, Phys. Rev. D 53, 1378 (1996); 
R. Akhoury and I. Rothstein, Phys. Rev. D 
54:2349 (1996); A. Leibovich and I. Roth- 
stein, Phys. Rev. D 61:074006 (2000); A. Lei- 
bovich, I. Low and I. Rothstein, Phys. Rev. 
D 61:053006 (2000); 

12. G. Bauer, S. Fleming and M. Luke, UTPT- 
00-03, hep/ph0005275. 

13. G. Balzereit, W. Kihan and T. Mannel, Phys. 
Rev. D 58:114029 (1998). 

14. A. Andrasi and J. G. Taylor, Nucl. Phys. B 
350:73-81 (1991). 

15. U. Aglietti, work in progress. 



Figure captions 

Fig.l (above): Integration region for the soft 
gluon in the plane (e, l±) with the Ag regulariza- 
tion. We have taken k+ = 0.3 GeV and As = 1 
GeV. 

Fig. 2 (below): Integration region for the soft 
gluon in the plane (e, l±) with the A^ rcgulariza- 
tion. We have taken fc+ = 0.3 GeV and A^ = 1 
GeV. The domain Dg is to the left of the dashed 
line, AD is to the right. 
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